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ABSTRACT 

Based on 19 high-resolution N-body/gasdynamical galaxy formation simulations in the ACDM cosmology 
it is shown, that for a galaxy like the Milky Way, in addition to the baryonic mass of the galaxy itself, about 
70% extra baryonic mass should reside around the galaxy (inside of the virial radius), chiefly in the form of hot 
gas. Averaging over the entire field galaxy population, this "external" component amounts to 64-85% of the 
baryonic mass of the population itself. These results are supported by the recent detection of very extended, 
soft X-ray emission from the halo of the quiescent, massive disk galaxy NGC 5746. Some of the hot gas may, 
by thermal instability, have condensed into mainly pressure supported, warm clouds, similar to the Galactic 
High Velocity Clouds (HVCs). Based on an ultra-high resolution cosmological test simulation of a Milky 
Way like galaxy (with a gas particle mass and gravity softening length of only 7.6xl0 3 IT 1 Mq and 83/r'pc, 
respectively), it is argued, that the hot gas phase dominates over the warm gas phase, in the halo. Finally, an 
origin of HVCs as "leftovers" from filamentary, "cold" accretion events, mainly occurring early in the history 
of galaxies, is proposed. 

Subject headings: cosmology: theory - galaxies: formation - methods: numerical 
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1. INTRODUCTION 

In the ACDM cosmology a disk galaxy like the Milky Way 
will have a virial mass of ~8xl0 n M & . For a universal baryon 
fraction ~ 0. 15, one would expect the baryonic mass of the 
Milky Way to be ~ 1 .2x 1 1 1 M© , assuming that the bulk of the 
baryonic material inside of r vlI (~ 250 kpc) is deposited onto 
the central galaxy. However, the baryonic mass (stars + cold 
gas) of the Milky Way is found to be just — 6xl0 10 M Q (e.g., 
Dehnen & Binney 1998; Sommer-Larsen & Dolgov 2001), so 
an amount of baryonic mass, as large as the actual mass of the 
Milky Way itself, appears to be "missing" (e.g., Silk 2004, 
Mailer & Bullock 2004). 

On the other hand, on theoretical grounds it has been known 
for long that galaxies (not only in groups and clusters, but also 
in the "field") should be embedded in extended haloes of hot 
gas (e.g., White & Frenk 1991; Sommer-Larsen 1991), and 
evidence for hot, dilute gas in the Galactic halo is quite strong 
(e.g., Sembach et al. 2003). Searches for X-ray emission from 
the haloes of external, quiescent disk galaxies have until very 
recently proved unsuccessful (as opposed to star-burst galax- 
ies, e.g., Strickland et al. 2002), which has been taken as an in- 
dication that the baryonic mass of such hot haloes is insignifi- 
cant, due to, e.g., strong AGN driven hot gas outflows at some 
point during galaxy formation (e.g., Benson et al. 2000). 

Based on recent cosmological simulations of disk galaxy 
formation (not invoking violent AGN feedback) Toft et al. 
(2002) showed, however, that the X-ray null detections were 
to be expected. Moreover, they showed that the X-ray lumi- 
nosity of disk galaxy haloes is expected to be a very steep 
function of the characteristic circular speed, roughly as Lx oc 
VJ. Very recently Pedersen et al. (2006) observed the mas- 
sive, quiescent, isolated and edge-on, disk galaxy NGC5746 
(V c = 305 ± 7 km/s), and detected hot halo soft (kT -0.4 keV) 
X-ray emission at the level predicted by the numerical models 
(see also Rasmussen et al. 2006 for more detail). As the mass 
of hot gas inside of r wa for such a galaxy is predicted to be 
— 80% of that of the central galaxy, and the total "external" 
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baryonic mass fraction — 1 10% (section 3), it is clearly of in- 
terest to estimate the global external baryonic mass-fraction, 
averaged over the entire field galaxy population. This is the 
aim of this Letter: in sec. 2 the code and simulations are 
briefly described, the results obtained are presented in sec. 3, 
and discussed in sec.4. Issues like formation of warm, mainly 
pressure supported clouds ("High Velocity Clouds") through 
thermal instability of hot halo gas, are also addressed 1 . 

2. THE CODE AND SIMULATIONS 

The code used for the simulations is a significantly im- 
proved version of the TreeSPH code, which has been used pre- 
viously for galaxy formation simulations (Sommer-Larsen, 
Gotz & Portinari 2003, SLGP03). The main improvements 
over the previous version are: (1) The "conservative" entropy 
equation solving scheme suggested by Springel & Hernquist 
(2002) has been adopted. (2) Non-instantaneous gas recy- 
cling and chemical evolution, tracing 10 elements (H, He, C, 
N, O, Mg, Si, S, Ca and Fe), has been incorporated in the 
code following Lia et al. (2002a,b); the algorithm includes 
supernovae of type II and type la, and mass loss from stars 
of all masses. (3) Atomic radiative cooling depending both 
on the metal abundance of the gas and on the meta-galactic 
UV field, modeled after Haardt & Madau (1996) is invoked, 
as well as simplified treatment of radiative transfer, switch- 
ing off the UV field where the gas becomes optically thick to 
Lyman limit photons on scales of — 1 kpc. 

For the present project, the formation and evolution of 15 
individual galaxies, known from previous work to become 
disk galaxies at z=0, was simulated with the above, signifi- 
cantly improved, TreeSPH code. At least two different numer- 
ical resolutions were used to simulate each galaxy. Moreover, 
many of the galaxies were also simulated with different phys- 
ical prescriptions for the early (z > 4) star-bursts (and related 
SNII driven energy feedback) found previously to be required, 
in order to produce realistic disk galaxies. The 15 galaxies 
were selected from a cosmological LCDM simulation (see be- 
low) to represent "field" galaxies, by requiring each galaxy to 

1 For images, and an HVC animation, see http://www.tac.dk/jslarsen/HVC 
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FIG. 1 . — For 19 high-resolution disk galaxy simulations are shown: total 
external baryon mass-fractions (solid squares), and contributions from hot 
gas (solid circles), satellite galaxies + outer halo stars (solid triangles), and 
inner halo stars (open circles) ; linear fits to the various components are shown 
by dotted lines. Also shown, by open symbols, are results for two elliptical 
galaxies (inner halo stars included in the central galaxy, hence not shown). 

be at least 1 Mpc from any galaxy group, and at least 0.5 Mpc 
away from any larger galaxy, at z=0. The disk galaxies span a 
range of characteristic circular speeds of V c ~ 100- 330 km/s, 
and a range of virial masses of 6xl0 10 to 3xl0 12 M Q . 

The galaxies (galaxy DM haloes) were selected from a cos- 
mological, DM-only simulation of box-length 10 /i _1 Mpc (co- 
moving), and starting redshift z,=39. The adopted cosmology 
was the flat ACDM model, with (Q M , ft A )=(0.3,0.7). 

Mass and force resolution was increased in Lagrangian re- 
gions enclosing the galaxies, and in these regions all DM 
particles were split into a DM particle and a gas (SPH) par- 
ticle according to an adopted universal baryon fraction of 
/i,=0.15, in line with recent estimates. In this paper, only 
results of high-resolution simulations, consisting of at least 
1.5xl0 5 SPH+DM particles, will be presented (with typi- 
cal numbers in the range 2-3xl0 5 ). Comparison to sim- 
ulations of lower and very high resolution, respectively, 
will be discussed, though. For galaxies of V c < 150 km/s, 
'Wgas='w*= 9.1xl0 4 and wdm=5.2x10 5 /i _1 Mg. Moreover, 
gravitational (spline) softening lengths of e gas =e*=190 and 
£dm=340 /r'pc, respectively, were adopted. For galaxies of 
V c > 150 km/s, m gas =OT st = 7.3xl0 5 and m DM =4.2xl0 6 /r'M©, 
and e gas =e*=380 and 6dm=680 /r'pc. In addition, two galax- 
ies of V c =180 and 244 km/s, respectively, were re-simulated 
with the smaller particle masses and gravity softenings above. 
These very high resolution runs consisted of 1.2 and 2.2 mil- 
lion particles. The gravity softening lengths were fixed in 
physical coordinates from z=6 to z=Q, and in comoving co- 
ordinates at earlier times. 

A Kroupa IMF was used in the simulations, and early rapid 
and self -propagating star-formation (sometimes dubbed "pos- 
itive feedback") was invoked (SLGP03). Finally, in order 
to enable some reuse of previous work, values of h=0.65 
and (78=1.0 were employed in the cosmological simulations. 
These values are slightly different from the h ~ 0.7 and <j% ~ 
0.9, currently favored. To check the effect of this, one simu- 
lation was undertaken with h = 0.7 and a% = 0.9. 

3. RESULTS 
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FIG. 2. — External fractions for 3 galaxies, simulated at high resolution 
with various star-formation and feedback prescriptions, versus total baryonic 
mass fractions within the virial radius (symbols as in Fig. 1 ; results for the 
same galaxy are connected by lines). For a given galaxy, as SNII feedback 
weakens, /j, v j r systematically increases. 

A detailed analysis of the simulation results, in general, 
will be given elsewhere, e.g., Portinari & Sommer-Larsen 
(2006); here emphasis will be just on "external" baryonic 
mass-fractions: For each disk galaxy simulation the distri- 
bution of baryonic mass, at z=0 and inside of the virial ra- 
dius, was classified as follows: 1) disk+bulge stars (R < 
20(V c /220km/s) kpc, |z| < 5(V e /220km/s)kpc), 2) inner halo 
stars (r < 20(V c /220km/s) kpc; not in the above region), 3) 
stellar satellites/outer halo (r < r^, and not in the above re- 
gions), 4) gas in the galaxy (mainly "cold": T < 3xl0 4 K; 
same region as for the disk+bulge stars), 5) outer hot gas 
(T > 3xl0 4 K; r < r v ; r , outside disk+bulge region), and 6) 
outer cold gas (r < r wlI ; outside disk+bulge region). The hot 
gas is typically at T ~ r v ; r , and hence, in general, much hotter 
than 3xl0 4 K. 

The simple estimate of the galactic baryonic mass would be 
the combined mass of components 1 and 4. In Fig.l is shown, 
for 19 high-resolution disk galaxy simulations, the mass of the 
other components relative to this baryonic galaxy mass. Com- 
ponents 3 and 6 have been combined, since most of the "exter- 
nal" cold gas is located in satellite galaxies. As can be seen, 
the importance of the three resulting "external" components 
increases with V c : at V c ~100 km/s, the external baryonic 
mass-fraction is only about 20%, at V c ~220 km/s (like the 
Milky Way) this fraction increases to ^70%, and at V c ~300 
km/s it is about a factor 1.1. At all V c , the hot gas is the domi- 
nant external component. Before inferring the implications of 
these results, their dependence on the physics invoked in the 
simulations, as well as on numerical resolution, should be as- 
sessed. To this end is compared in Fig. 2 the results of various 
simulations of three galaxies, all run at high resolution, with 
different prescriptions for SNII feedback. For a small galaxy, 
of V c ~ 1 10 km/s, is shown results for four different simula- 
tions: two with early, self -propagating star-formation (SPSF), 
and associated strong feedback (for two different threshold 
densities «H.e.iow=0. 1 and 0.25 cm" 3 ; SLGP03), one with early 
fast star-formation, but no SPSF (this reduces the early feed- 
back considerably), and one with no early fast star-formation 
at all (this corresponds to the very weak feedback case). Fig. 2 
shows that a) the stronger the feedback the smaller is the re- 
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FIG. 3. — Galaxy baryonic mass (disk+bulge+cold gas), (solid squares), 
total virial mass (solid triangles), and hot gas mass (solid circles) versus Vc. 
Dotted lines show power law fits. 

suiting baryon fraction inside of the virial radius, and b) the 
quantities shown in Fig.l are quite robust to such a dramatic 
change of feedback physics (the two lower feedback simula- 
tions result in galaxies with smaller specific angular momenta 
and stellar disks, which are kinematically too hot compared 
to observed disks — hence the results for these have not been 
included in Fig.l). Also shown are data for SPSF simula- 
tions, with the above two thresholds, for V c ~205 km/s and 
a V c ~305 km/s galaxies, respectively. As can be seen, these 
results are quite robust as well. 

To test for resolution effects, results for 7 galaxies with V c in 
the range 170-250 km/s, run at low resolution (the main one of 
SLGP03) and high resolution, but identical feedback descrip- 
tions, were compared. It was found that the hot gas fraction 
decreased by 28 ±4%, and the total "external" baryonic mass 
fraction by 18±6%, increasing the mass resolution by a fac- 
tor of eight, and the force resolution by a factor of two. For 
a single galaxy of V c ~305, the corresponding numbers were 
23 and 13%. To test the effects of going to very high resolu- 
tion, results for two galaxies, of V e =180 and 244 km/s, were 
compared at z ~ 1 and 2, respectively (results at z=0 were not 
available, but at both redshifts gas cool-out and star formation 
is very well underway). It was found that in going from high 
to very high resolution, the hot gas fraction is reduced only by 
about 5% more, and the total external baryon fraction did not 
change at all. On basis of the above, only simulations run at 
high resolution were used in the present analysis. 

4. DISCUSSION 

To estimate the global importance of the external baryons 
in the mass budget of disk galaxies, the global distribution 
of disk galaxy circular velocities (the "velocity function") is 
required. Galaxy velocity functions have been estimated by, 
e.g., Gonzalez et al. (2000) and Desai et al. (2004). Most 
published velocity functions refer to the general morphologi- 
cal mix of disk, lenticular and elliptical galaxies, in the field 
or in clusters, but Gonzalez et al. (2000) give an expression 
for field disk galaxies, viz. 



*d(V c ) dV c = 9 d . ex P 
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where = 2.0 ±0.4 10" 2 h 3 Mpc" 3 , 0= 1.3 ±0.18, n=2.5 
and V c * = 247 ± 7 km/s. Next, one has to determine the re- 



lation between V c and the combined baryonic mass of the 
disk+bulge stars and gas in the galaxy. In Fig. 3 this relation 
is shown for the 19 high resolution simulations. The relation 
is very well fitted by a power law 



M oa i = (5.0±0.1)xl0 1 
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(3.00±0.07) 
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Thirdly, the external baryonic mass-fraction, k = 
(A^bar.extemai/^gai), is parameterized as a function of V c 
by a linear fit to the data shown in Fig.l, resulting in 



k(V c ) = (1.00 ±0.07) 



220 km/s 



-(0.31 ±0.02) . (3) 



The fit is valid in the range V c ^100-330 km/s, and extrapo- 
lates to zero at V c =68 km/s; in the following n=0 is assumed 
for V c < 68 km/s, to obtain a lower limit to the global bary- 
onic external mass-fraction (actually, at V c <100 km/s, k is 
expected to start increasing with decreasing V c , as SNII feed- 
back increasingly suppresses gas cool-out and star-formation: 
for a ^3xl0 5 particle simulation of a V c =41 km/s galaxy, an 
external baryon fraction of ~ 1.4 was found). 

The globally averaged external baryonic mass-fraction can 
now be expressed as 



K = ■ 



Jo ' mm KjVc) M gal (V c ) -$(V C ) dV c 
Jo ycmax M gal (V c ) ¥(V C ) dV c 



(4) 



where Vc iinax is the maximum circular speed of disk galaxies. 
Using eqs.[l]-[3], and assuming V cmax =350 km/s (disk galax- 
ies of such circular speeds are certainly observed), k=0.64 is 
obtained. Increasing V c -. max to 500 km/s (which must be con- 
sidered an observational upper limit), results in an increase 
of R to 0.72. Fitting to the various external components indi- 
vidually, shows that hot gas contributes 70% to R, inner halo 
stars 5%, and outer halo stars, satellite stars and cold gas in 
satellites 25%. Assuming that the above results on external 
mass fractions also apply to field elliptical/lenticular galaxies 
(see Fig.l), and applying the field galaxy velocity function of 
Desai et al. (2004), increases the estimate of R to 0.72 and 
0.85, for V c , max =350 and 500 km/s, respectively (in principle, 
very strong AGN feedback, not included in the current simu- 
lations, could lower hot gas fractions of E/S0 galaxies; on the 
other hand E/SOs are (contrary to disk galaxies), in general, 
observed to be embedded in hot gas haloes). Finally, since, in 
particular, the larger galaxies host a number of smaller galax- 
ies within r y -„, the velocity functions should actually be re- 
duced accordingly at the low V c end. Hence, the above R esti- 
mates are, in fact, lower limits. 

The above results, taken at face value, provide the explana- 
tion why only about half of the baryonic mass, inside of the 
virial radius of galaxies like the Milky Way, is actually found 
in stars and cold gas in the galaxies. Given the ramifications 
of this result, it is important to test its robustness: Although 
fb ~0.15 is indicated by standard big bang nucleosynthesis, 
various CMB results and direct measurements of ft, in galaxy 
clusters, lower values are still not excluded. To test the effect 
of adopting ft, =0. 10, five galaxies were re-simulated with this 
value, at the resolution of SLGP03. The results of these 5 
simulations were compared to a sample of 39 simulations run 
with f h =0.15 at the resolution of SLGP03. It is found that 
hot gas fractions, as well as total external baryonic mass frac- 
tions, increase by 30-50% (at a given V c ) going from fb =0.15 
to 0.10. This is mainly due to the decreased efficiency of hot 
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FIG. 4. — Mass distribution of warm clouds at r <50 kpc for the ultra-high 
resolution simulation at f=0.3 Gyr (dashed line), and 0.7 Gyr (solid line), as 
well as for the very high resolution simulation at t=0J Gyr (dotted line). 

gas cool-out, and the resulting increase of the virial radius at 
a given circular speed of the central disk galaxy. 

To test the effect of adopting h=0.7 and <ro=0.9, rather than 
the 0.65 and 1.0 used in the main simulations, one simula- 
tion was run with the former parameters at the resolution of 
SLGP03. Comparing to the above 39 simulations it was found 
that the hot gas, as well as, total external fraction increase by 
~5-10%. This is likely mainly due to the Hubble time being 
7% less, allowing less time for hot gas cool-out. 

4.1. Thermal instability of hot halo gas, and HVCs 

The hot, dilute halo gas may be susceptible to thermal insta- 
bility, causing the formation of a two-phase medium, consist- 
ing of a hot phase (T ~ 10 6 K) and a warm phase (T ~ 10 4 K), 
in approximate pressure equilibrium (e.g., Burkert & Lin 
2000; Mailer & Bullock 2004, Kaufmann et al. 2005). This 
may not change the general picture, since warm clouds, once 
formed, may quickly be destroyed again due to various phys- 
ical processes. It is, however, important to estimate the frac- 
tion of halo gas in warm clouds, at any given time. To es- 
timate this for a realistic, cosmological simulation very high 
numerical resolution is required, and the following approach 
was adopted: for a "standard" high-resolution simulation of 
a Milky Way like galaxy (V c =224 km/s, M vir =8xlO n M , 
r v i r =256 kpc), at z=0.4, each gas particle, at r <300 kpc, was 
split into 12 equal mass gas particles, to achieve very high 
resolution of the gas phase inside of r v i r . The simulation was 
then continued for 0.2 Gyr, whence each gas particle, at r <50 
kpc, was again split into 8 equal mass gas particles, to achieve 
ultra-high gas phase resolution in the region where most of the 
warm clouds reside (see below). The resulting ultra-high reso- 
lution gas particle mass is 7.6xl0 3 h~ l M Q , and gravity soften- 
ing length eg as =83/! _1 pc. The simulation was then continued 
for an additional 0.5 Gyr. In Fig.4 is shown, for r < 50 kpc, 
the distribution of warm cloud masses, at f=0.3 and 0.7 Gyr. 
The integrated warm cloud mass decreases somewhat with 
time (-20%), to ~5xlO 7 M at t =0.7 Gyr, or about 2% of the 
total mass of hot halo gas at r < 50 kpc. Also shown, is the 



mass distribution of warm clouds, at r < 50 kpc, in the very 
high resolution simulation, at f=0.7 Gyr. To the cloud resolu- 
tion limit of this simulation, M c i res ~2xlO 6 M , the agreement 
between the two simulations is quite good. Moreover, with a 
resolution limit of M c i , es ^3xlO 5 M for the ultra-high reso- 
lution simulation, the peak at M c i ~7-8xlO 5 M appears well 
resolved. It hence seems unlikely that the mass contribution 
from warm clouds of mass less than about 3xlO 5 M is sig- 
nificant; such clouds will in any case quickly be destroyed by 
various physical processes, e.g., Mailer & Bullock (2004). 

From the very high resolution simulation, the mass in warm 
clouds at 50< r <256 kpc, is about 75% of that within 50 kpc. 
Assuming this ratio, the total mass of warm clouds, at r < r v ; r , 
down to M c i ~3xlO 5 M , estimates to ~1O 8 M . It seems 
unlikely that inclusion of warm clouds of M c i <3xlO 5 M 
(which can not be resolved by the current simulation) will 
increase this estimate by more than a factor of a few. Put- 
man (2006) estimates the total gas mass of the HVC system 
of the Milky Way to ^4-6x1 8 M , assuming that the HVCs 
are distributed to distances of ^60 kpc (a reasonable assump- 
tion according to the present findings). This is in fair agree- 
ment with, though somewhat larger than, the above estimates. 
Moreover, the hot gas in the halo of the simulated Milky Way 
like galaxy, has n H ~10" 3 5 ,10 -4 and 10" 4 5 cm" 3 at r ~10, 50 
and 100 kpc. This is in reasonable agreement with the lower 
limit of ~ 1 -4 cm" 3 , deduced by Quilis & Moore (200 1 ) from 
the observed head-tail structure of many Galactic HVCs (also 
displayed by many of the simulated warm clouds) 1 . 

The warm clouds are mainly confined by the pressure of 
the ambient halo gas, though some of the most massive dis- 
play high density cores, and have a ratio of gravitational to 
thermal energy of about 0.5. Moreover, they typically appear 
to be seeded by warm filamentary structures, which are left- 
overs from earlier "cold" accretion (e.g., Birnboim & Dekel 
2003; Keres et al. 2005; Sommer-Larsen 2005 — details will 
be given in Putman & Sommer-Larsen 2006). No additional 
warm clouds appear to form, despite that the necessary con- 
dition for onset of thermal instability, viz. t\ < r coo i, is satis- 
fied everywhere in the hot halo gas (t> is the sound crossing 
time, which is taken to be ~ 2/isph/c s , where /isph is the local 
SPH smoothing length and c s is the sound speed: r coo i = | 
is the timescale for radiative cooling). The mass of hot gas 
within r v j r is ~2.3xlO lo M , so the ratio of warm cloud mass 
to hot gas mass is at most a few percent. Although fully cos- 
mological simulations at even higher resolution are required 
to address these issues completely, the above results strongly 
suggest that the hot gas phase is the dominant one in the halo. 
It is also worth noting that Pedersen et al. (2006) found that 
the predicted halo X-ray luminosity and surface brightness 
profile of large disk galaxies, simulated at a resolution which 
was, though high, not sufficiently high to resolve warm halo 
clouds, matches the observed X-ray properties of NGC5746 
very well. This would not be expected if a substantial fraction 
of the halo gas was in the warm phase. 

I have benefited from discussions with J. Binney, A. Dekel, 
R. Bower, J. Fynbo, L. Hernquist, T. Kaufmann, A. Loeb, 
L. Portinari, M. Putman and J. Silk. 

The Dark Cosmology Centre is funded by the DNRF. 
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